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ABSTRACT 


Representation and identification of a nonlinear aeroelastic pitch-plunge system as 
a model of the Nonlinear AutoRegressive, Moving Average exogenous (NARMAX) class 
is considered. A nonlinear difference equation describing this aircraft model is derived 
theoretically and shown to be of the NARMAX form. Identification methods for NARMAX 
models are applied to aeroelastic dynamics and its properties demonstrated via continuous- 
time simulations of experimental conditions. Simulation results show that (1) the outputs 
of the NARMAX model closely match those generated using continuous-time methods, 
and (2) NARMAX identification methods applied to aeroelastic dynamics provide accurate 
discrete-time parameter estimates. Application of NARMAX identification to experimental 
pitch-plunge dynamics data gives a high percent fit for cross-validated data. 
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maximum error lag 
maximum input lag 
maximum output lag 
record length 
validation record length 

Nonlinear AutoRegressive, Moving Average exogenous 

prediction error identification 

signal-to-noise ratio, dB 

standard deviation 

sample time 

controlled or exogenous input 
free-stream velocity 
weighted least-squares 

nondimensional distance between elastic axis and center of mass 

output of the continuous time simulation 

simulated output 

measured output 

percent quality of fit 

nonlinear mapping 

pitch angle 

pitch velocity 

pitch acceleration 

control surface deflection 

Kronecker impulse function 

density of air 


INTRODUCTION 

System identification, or mathematical modeling, is the process of developing or 
improving a mathematical representation of a physical system based on observed data. 
System identification is a critical step in aircraft development, analysis, and validation for 
flight worthiness. 

One such application of system identification in the flight-test community is for the 
analysis of aeroelasticity. The analysis of aeroelasticity is concerned with the interaction 
of inertial, structural, and aerodynamic forces (ref. 1). Previous approaches have 
modeled aeroelasticity with linear time-invariant (LTI) models. These linear models have 


been successful in providing approximate estimates of an aircraft’s response to gust, 
turbulence, and external excitations. However, when aircraft speeds increase to high 
subsonic ortransonic Mach numbers, linear models no longer provide accurate predictions 
of the aircraft’s behavior. Some of the behavior that cannot be modeled linearly includes 
transonic dip, airflow separation, and shock oscillations, which can induce nonlinear 
phenomena such as limit cycle oscillations (LCO) (refs. 2 and 3). The onset of LCOs has 
been observed on several aircraft, such as the F-16C (Lockheed Martin Corporation, 
Bethesda, Maryland) or F/A-18 (McDonnell Douglas Corporation, St. Louis, Missouri, 
now The Boeing Company, Chicago, Illinois), and cannot be modeled properly as an LTI 
system (ref. 4). This limitation has necessitated the application of nonlinear identification 
techniques to accurately model LCO dynamics. 

Over the past several decades, significant achievements have been made in several 
areas of nonparametric nonlinear system identification (e.g., refs. 5-7). Recent work in the 
aerospace community has attempted to address these nonlinear aeroelastic phenomena 
using Volterra kernel methods (ref. 8). These methods provide a convenient means of 
characterizing LCOs but suffer from a highly over-parameterized model description and 
do not lend themselves to efficient control synthesis. 

Parametric representations of nonlinear systems typically contain a small number 
of coefficients that can be varied to alter the behavior of the equation and may be linked 
to the underlying system. Leontaritis and Billings (refs. 9 and 10) have proposed the 
Nonlinear AutoRegressive, Moving Average exogenous (NARMAX) structure as a general 
parametric form for modeling nonlinear systems. NARMAX models describe nonlinear 
systems in terms of linear-in-the-parameters difference equations relating the current 
output to (possibly nonlinear) combinations of inputs and past outputs. It is suitable for 
modeling both the stochastic and deterministic components of a system and is capable 
of describing a wide variety of nonlinear systems (refs. 11 and 12). This formulation 
yields compact model descriptions that may be readily identified and may afford greater 
interpretability. NARMAX models have been successfully demonstrated for modeling the 
input-output behavior of many complex systems such as ones found in engineering and 
biology (refs. 13 and 14). 

Currently, development and test of aircraft takes many years and considerable 
expenditure of limited resources. One reason for lengthy development time and costs 
is that many models (and hence control strategies) need to be developed throughout 
the flight envelope. The power of parametric nonlinear identification techniques (i.e. , 
NARMAX models) is that they can describe complex aeroelastic behavior over a large 
operating range. Consequently, models are provided that can be more robust and reduce 
development time. 

Although the NARMAX structure is well-suited to modeling the input-output behavior 
of an aeroelastic system, this method has not been investigated by the flight-test 
community to date. Therefore, the objectives of this report are to (1 ) theoretically analyze 
a nonlinear pitch-plunge model of aircraft dynamics to derive its NARMAX representation, 
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(2) assess the applicability of this nonlinear model for the identification of aerospace 
systems, and (3) investigate the suitability of NARMAX identification methods applied to 
aircraft dynamics. 

The results show that the NARMAX model class provides an ideal framework for 
modeling the input-output behavior of a nonlinear pitch-plunge model describing aircraft 
dynamics. Identification results illustrate that methods for identification of NARMAX 
models are well-suited for identifying aircraft dynamics. Analysis of experimental data 
using NARMAX identification techniques provides a parameter set that explains the 
input-output data well. Overall, this report contributes to the understanding of the use of 
parametric identification techniques for modeling of aerospace systems. 

The organization of this report is as follows. The NARMAX model structure is described 
in the section titled “NARMAX Model.” In the section titled “Nonlinear Pitch-Plunge Model 
of Aircraft Dynamics,” a continuous-time representation of a nonlinear pitch-plunge model 
describing aircraft dynamics is given, while its NARMAX representation is derived in the 
section titled “Theoretical Analysis.” The section titled “Validation of NARMAX Pitch- 
Plunge Model” illustrates the results of simulating this NARMAX representation of pitch- 
plunge dynamics. In the section titled “NARMAX Identification of Pitch-Plunge Model,” 
the applicability of NARMAX identification to this model representation via simulations 
of experimental condition is assessed. The section titled “Identification of Experimental 
Pitch-Plunge Data” presents the results of identifying experimental wind tunnel data, and 
the “Discussion” section provides a discussion of the major finds. Lastly, conclusions and 
significance of the results are given in the “Conclusion” section. 

NARMAX MODEL 


The NARMAX structure is a general parametric form for modeling nonlinear systems 
(ref. 9). This structure describes both the stochastic and deterministic components of 
nonlinear systems. Many nonlinear systems are a special case of the general NARMAX 
structure (ref. 12). The NARMAX structure models the input-output relationship as a 
nonlinear difference equation of the form shown in equation (1 ): 


?(n) = / / [z(n-1),---,z(n-n z ),u(n),---,u(n-n u ),e(n-1),---,e(n-n e )] + e(n) (1) 


/denotes a nonlinear mapping; u is the controlled or exogenous input; z is the measured 

output; and e is the uncontrolled input or innovation. The nonlinear mapping, f 1 , can be 
described by a wide variety of nonlinear functions such as a tanh(-) or splines (i.e. , hard 
nonlinearities) (refs. 11 and 12). For simplicity, nonlinearities are only considered that can 
be described by a polynomial expansion. This class of nonlinear difference equations 
describes the dynamic behavior of a system as a linear and/or nonlinear expansion of 
the input, output, and error. These equations may include a variety of nonlinear terms, 

such as terms raised to an integer power (e.g., u 3 (n- 5)); products of present and 
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past inputs (e.g., u(n)u(n-3)Y, past outputs (e.g., z 2 (n-2)z(n-7)); or cross-terms 

(e.g., u{n-3)z 2 {n-4)). This system description encompasses many forms of nonlinear 
difference equations that are linear-in-the-parameters. Thus, there are no problems with 
local minima. 

NONLINEAR PITCH-PLUNGE MODEL OF AIRCRAFT DYNAMICS 

O’Neil et al. (refs. 15 and 16) developed a nonlinear pitch-plunge model (fig. 1) 
describing aircraft aeroelastic dynamics. Figure 1 characterizes aeroelastic wing dynamics 
for experiments performed on the Texas A&M University (College Station, Texas) testbed 
(ref. 16). This model provides a relationship between control surface deflection as input, 
and pitch-plunge displacement and velocity as outputs, of a single-input multiple-output 
nonlinearity followed by a simple integrator. 

The model presented in figure 1 has been derived from the governing equations of 
motion for aeroelastic systems, as shown in equation (2): 
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where h denotes plunge motion; a is the pitch angle; x a is the nondimensional distance 
between elastic axis and center of mass; m is the wing mass; l a is the mass moment 
of inertia of the wing about the elastic axis; b is the semichord of the wing; c h and c a 
are the plunge and pitch structural damping coefficients, respectively; k h and k a are 
the plunge and pitch structural spring constants, respectively; and L and M are the 
aerodynamic lift and moment, respectively. 
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Figure 1. System structure assumed for modeling and identification of pitch-plunge 
aeroelastic dynamics. 


Typically, quasi-steady aerodynamic forces and moments are assumed which can 
be modeled, as shown in equation (3): 
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where p denotes density of air; U is the free-stream velocity; c m and c, are the 


'L 


a a 

moment and lift coefficients per angle of attack, respectively; c m and c y are the 

moment and lift coefficients per control surface deflection, respectively; {3 is the control 
surface deflection; and a is the nondimensional distance from midchord to elastic axis. 


Although several classes of nonlinear mappings for stiffness contributions, k a (oc), 

have been investigated for open loop dynamics of aeroelastic systems (refs. 17-20), the 
work of O’Neil etal. (refs. 15 and 16) demonstrated that a polynomial mapping of the form 
shown in equation (4) describes the behavior of this testbed well: 
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Equations of motion are derived by combining equations (2) and (3) to yield 
equation (5): 
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The model form presented in figure 1 is derived by transforming equation (5) to give 
equation (6): 


x = yx) + g(x)^ (6) 

where x = [x 1 x 2 x 3 x 4 ] r = [/? a h a] 7 "; p = U 2 \ and f and g(x) are as given in 
equation (7): 
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The supplementary variables k, , i = 1 , 2, 3, 4, and g y - , / = 3, 4, are provided in table 1 in 
relationship to the aeroelastic parameters given in equations (2) and (3). 
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Table 1 . Supplementary system variables. 
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Transformation of equation (5) to equation (6) and introduction of supplementary 
variables (table 1) provides the simple model description presented in figure 1. The 
nonlinear mappings for this model description are given as equations (8) and (9): 


fg() = - k 2 h ( n ) 
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Note that this system (fig. 1) can be described in terms of pitch-plunge displacement or 
velocity. Here, pitch-plunge is chosen in terms of velocity because (1) it offers a model 
description with lower dynamic order, and (2) velocity feedback models are often used for 
vibration suppression. 


Discrete-Domain Approximations 

Many methods exist for discretization of continuous-time systems and signals. Most 
commonly used are Euler’s forward, Euler’s backward, or Tustin’s method (also known 
as the bilinear transformation method) (refs. 21 and 22). Each has its advantages as 
well as disadvantages. While Tustin’s and Euler’s backward methods provide a superior 
approximation to a continuous-time signal, the goal in system identification is not (directly) 
signal reproduction but model estimation. 

Tustin’s method provides excellent signal estimation, but its use for modeling a pure 
integrator yields system descriptions that can be overly complex (e.g., redundant terms). 
Euler’s backward method also provides good signal estimation but provides a model that 
is not intuitive. Models based on this approximation would include the current output as 
one of the model terms, leading to an algebraic loop. Although Euler’s forward method 
is well-known to be unstable, this is only true if the sampling rate is not sufficiently large. 
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For identification purposes this does not pose a concern since the signals need to be 
sampled at least twice Nyquist and, hence, stability is achieved. Generally, the rule of 
thumb is to sample a signal at least 4 to 10 times the highest known (or suspected) 
system dynamics (ref. 23). For the pitch-plunge system under investigation, Euler’s 
forward method provides a model description that is both stable and intuitive. Moreover, 
all three methods converge to similar accuracies for sufficiently large sampling rates. For 
these reasons Euler’s forward method is chosen to model the system dynamics. 


The pitch-plunge model is given in continuous-time. This section shows how the 
model can be converted to discrete-time and rewritten as a NARMAX model. To do so, 
note that the two nonlinearities can be decoupled and analyzed separately since they 
yield two separate model descriptions for pitch and plunge velocity. 

Euler’s forward (explicit) method (ref. 21), as shown in equation (10) 


where T is the sample time, was used to approximate the continuous-time integrator, 
where x is replaced by a and h for pitch and plunge, respectively. 

The nonlinearities used for this analysis, to derive an input-output model of pitch 
and plunge, were given in equations (8) and (9). In addition, the models are assumed 
corrupted by output additive (measurement) noise, as shown in equation (11): 


where d(n) and h(n ) are the noise-corrupted outputs, a n f(n) and h n f(n) are the 
unmeasured noise-free outputs, and e^(/i) and e^(n) are the measurement noise. 

After collecting and combining terms, the overall nonlinear models were represented 
as nonlinear difference equations with 10 terms each, as shown in equation (12): 


THEORETICAL ANALYSIS 
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( 12 ) 


a{n) = yyu{n - 1) + y 2 h(n - 1) + y 3 a{n - 1) + y 4 a{n - 1) 2 + y 5 a{n - 1) 3 

+y 6 oc{n - 1) 4 + y 7 a{n - 1) 5 + y Q h(n - 1) + y Q u(n - 1) ■ + y 1 Q e d (n - 1) ■ + e & (n) 
h(n) = O^h(r) - 1) + 0 2 h{n - 1) + 0 3 a{n - 1) + 0 4 a{n - 1) 2 + 0 5 a{n - 1) 3 

+0 6 a{n - 1) 4 + d 7 a(n - 1) 5 + 0 Q a{n - 1) + 0 Q u{n - 1) + 0 1 {n - 1) + e /? (n) 

These are NARMAX models since they (1 ) include input-output terms that are combinations 
of linear and nonlinear integer powers and (2) are linear-in-the-parameters. Table 2 shows 
the relationship of discrete-time NARMAX parameters in equation (12) to the underlying 
continuous-time coefficients. 


Table 2. Theoretical relationship of NARMAX model parameter set to continuous-time 
system coefficients. 


NARMAX 

plunge 

coefficient 

Relationship to 
continuous- 
time coefficient 

NARMAX 

pitch 

coefficient 

Relationship to 
continuous- 
time coefficient 


1-737 

7i 

1 -Tt^ 

0 2 

Ta 1 

72 

Tb 1 

^3 

C\J 

i 2 

7 3 

Tb 2 

04 

CO 

£ 

7 4 

Tb 3 

05 

Ta 4 

7 5 

Tb 4 

06 

Ta 5 

7 6 

Tb 5 

07 

Ta 6 

7/ 

Tb 6 

08 

Ta 8 

7 8 

Tb, 

09 

S' 

CD 

7g 

CD 

0 io 

-(i-ra 7 ) 

7io 



11 


VALIDATION OF NARMAX PITCH-PLUNGE MODEL 


The accuracy of this system representation was validated by simulating the pitch- 
plunge model in continuous-time using Simulink® (The MathWorks, Natick, Massachusetts) 
(fig. 1). The nonlinearities used in this continuous-time simulation were the fifth-order 
power series described in equations (8) and (9). The parameters used in the simulation 
were typical values found in experiments and are given in table 3 (ref. 15). The system 
was excited using a 5-Hz chirp input. 


Table 3. Continuous-time system coefficients. 


CT coefficient 

Value 

b 

0.135 m 

span 

0.600 m 

k h 

2844.4 N/m 

c h 

27.43 N s/m 

P 

1 ,225 kg/m 3 

c i 

a 

6.28/rad 

c. 

'p 

3.358/rad 

° m a 

(0.5 + a)c, /rad 

a 

% 

-0.635/rad 


Output Accuracy 

To determine the validity of this NARMAX description model provided in equation (12), 
its response is simulated for a parameter set corresponding to those used for the 
continuous-time model. The input sequence was a 5-Hz chirp with a signal duration of 
30 s. The chirp input had an operating range between ±1 .0 rad, as shown in figure 2(a). 

Using a fifth-order nonlinearity, the frequency content of the signal at the output of 
the nonlinearity will be at least 25 Hz (5 times the 5-Hz chirp signal). To avoid internal 
aliasing, a sampling rate of 1 00 Hz ( T = 0.01 s) is selected, which is 4 times greater than 
the internal 25-Hz signal. 
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The simulated output, y{n), of the NARMAX description model was compared with 
the output of the continuous-time simulation, y(n), by computing the percent variance 
accounted for by the NARMAX model as the percent quality of fit (%QF) provided in 
equation (13): 


f 

%QF = 1 


v 




x 100 


where N is the record length. 


(13) 


Simulation Result 

The results of simulating the pitch-plunge model in continuous-time against the 
discrete-time NARMAX predictions are provided in figure 2. Figure 2 shows the simulation 
input (fig. 2(a)) and predicted velocity outputs of the NARMAX description models 
superimposed on top of the continuous-time outputs of the pitch (fig. 2(b)) and plunge 
(fig. 2(c)) velocity models. With over 99 %QF, the NARMAX outputs matched that of the 
continuous-time simulation with negligible error. 



Time, s 


080313 


(a). Control surface deflection. 


Figure 2. (a). Input to simulated pitch-plunge model in continuous-time and NARMAX 
description model; (b) and (c). Output of simulated pitch-plunge velocity model in 
continuous-time and NARMAX description model. Note that the two outputs are almost 
identical. 
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(b). Pitch velocity. 


(c). Plunge velocity. 


Figure 2. Concluded. 


NARMAX IDENTIFICATION OF PITCH-PLUNGE MODEL 

The utility of methods developed for identifying NARMAX models using sampled 
data from this continuous-time simulation were then assessed. An extended least-squares 
(ELS) algorithm (refs. 24-26) was used to identify model parameters. 

The NARMAX description of the pitch-plunge velocity models (eq. (12)) is described 
by past outputs that are linear-in-the-parameters. In the presence of output additive noise 
(eqs. (11) and (12)), these terms result in lagged values of disturbance terms that are 
also linear-in-the-parameters. If these lagged errors are not modeled, they induce a bias 
in the parameter estimates (refs. 22, 27, and 28). The ELS algorithm was implemented 
because it is designed to model lagged-error terms thereby providing unbiased parameter 
estimates. 

The fact that ELS may suffer from convergence problems (refs. 22, 29, and 30) is 
well-known. However, no prediction error identification (PEI) method is optimal. For the 
pitch-plunge velocity models, ELS is deemed the best estimation technique because it 
provides an unbiased estimate of model parameters (ref. 22). Other estimation techniques 
such as maximum likelihood (ML), instrumental variables (IV), weighted least-squares 
(WLS), etc., are difficult to implement and also have convergence problems (refs. 22, 30, 
and 31). For this reason ELS is chosen. 
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Analysis of NARMAX Model Parameters 

A Monte Carlo study of these NARMAX parameters (eq. (12) and table 2) was 
performed to assess their accuracy and variability using the ELS estimator. One thousand 
Monte Carlo simulations were generated in which the input-output realization was the 
same but had a unique Gaussian, white, zero-mean, noise sequence added to the output. 
The output additive noise amplitude was increased in increments of 5 dB, from 20 to 
0 dB signal-to-noise ratio (SNR). Parameter mean and standard deviation was computed 
from the 1000 estimates. The input used for this study was the same 5-Hz chirp signal 
previously described in the section titled “Output Accuracy.” 

For this study, the system order and structure were assumed to be known with the 
coefficient set in equation (12) and tables 2 and 3. The regressor matrix used by this 
algorithm was formed to contain only those columns (parameters) that corresponded to 
the theoretical analysis (eq. (12)). A reasonable assumption is made that the order and 
structure are known because the goal is to identify a model in the model class described 
by equation (12). 

In experimental settings, often only pitch-plunge displacement and/or acceleration 
signals are available. For a velocity model description, the velocity signal is required for 
identification. Therefore, pitch-plunge acceleration signals are numerically integrated to 
obtain a velocity profile. The estimation set consisted of N - 3000 data points sampled at 
T = 0.01 s. The estimated parameters were cross-validated with a fresh noise-corrupted 
output to compute the %QF of the predicted pitch and plunge velocity. The validation set 
consisted of N v = 3000 data points (refs. 22 and 32). 

Identification Results of Simulated Model 

Figure 3 shows the results of identifying this simulated model of pitch-plunge. The 
NARMAX parameters in this figure correspond to those given in table 2. This figure shows 
that the identified parameter values corresponded closely to those derived theoretically 
for all SNRs. Note that the mean value of parameters y 10 and 0 1O is not expected to be 

close to the theoretically computed value since they correspond to lagged-error terms. 
Lagged-error terms are difficult to identify accurately, even with high SNR, since the 
terms model the output additive noise, which is an unmeasurable stochastic process. 
This stochastic process is modeled (approximated) by a deterministic signal of prediction 
errors, which is only an (poor) estimate of the noise (refs. 22 and 28). 
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(a). Pitch parameters, mean and STD. 


(b). Plunge parameters, mean and STD. 


Figure 3. Monte Carlo study of pitch-plunge NARMAX model parameters. 5-Hz chirp input, 
Gaussian, white, zero-mean noise and A/= 3000. Ordinate: STD about mean. Abscissa: 
Output SNR = 20, 15, 10, 5 and 0 dB. (Note that the abscissa is shown in decreasing 
SNR, which corresponds to increasing noise intensity.) 


Figure 4 presents a result of cross-validation for a typical parameter set for this 
study. Figures 4(a) and 4(b) show a noise-corrupted output used for identification, and 
figures 4(c) and 4(d) show a predicted output superimposed on top of the noise-free 
output. The predicted output matched the measured output with over 98 %QF. 
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(c). Predicted pitch velocity: QF = 98.33%. (d). Predicted plunge velocity: QF = 98.54%. 


Figure 4. Cross-validation for typical identified NARMAX pitch (left) and plunge (right) 
models with N = 3000 and Gaussian, white, zero-mean output additive noise (0 dB 
SNR). Top panels: Measured outputs used for estimation. Bottom panels: Predicted 
output superimposed on top of noise-free output. 


IDENTIFICATION OF EXPERIMENTAL PITCH-PLUNGE DATA 

Lastly, the identification technique is assessed on experimental wing section data 
collected in the wind tunnel at the Department of Aerospace Engineering at Texas A&M 
University by the Aeroelasticity Research Group. The data analyzed for this study does 
not contain a flap control input but instead had an initial condition associated with plunge 
displacement. Data with a control input was unavailable for analysis. 

The velocity model descriptions given in equation (12) are in terms of a control input 
and zero initial conditions. For the present study the model was modified to reflect a lack 
of exogenous input and the presence of an initial condition. Equation (12) is reformulated 
as equation (14): 
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( 14 ) 


h(n ) = Oji(n - 1) + 0 2 h(n - 1) + 0 3 a{n - 1) + 0 4 a(n - 1) 2 + 0 5 a(n - 1) 3 

+0 6 a{n - 1) 4 + 0 7 a{n - 1) 5 + 0 Q a{n - 1) + 0 g 8(n) + 0 1o e^(n - 1) + e f 1 ( n ) 
a(n) = y^a(n - 1) + y 2 h(n - 1) + y 3 a(n - 1) + y 4 a(n - 1) 2 + y 5 a(n - 1) 3 
+ 7 6 a(n - 1) 4 + y 7 a{n - 1) 5 + y Q h{n - 1) + y Q e d (n - 1) + e d (n) 

where 8(n) is the Kronecker impulse function used to represent the onset of a plunge 
initial condition in discrete-time. Note that this model description can also be modified 
for use in analysis of data that contains both an initial condition and exogenous input by 
simply including a Kronecker impulse function in equation (12) or excluding both (time- 
series analysis) by removing the exogenous input term. 

Apparatus 

Data was collected on a unique wind tunnel test apparatus at the Department of 
Aerospace Engineering at Texas A&M University. This 2-ft by 3-ft closed-circuit low-speed 
wind tunnel allows a wing section to move in two degrees-of-freedom and can translate 
(plunge) and rotate (pitch). This apparatus allows the study of classical bending-torsion 
flutter. Structural response of the system is governed by springs attached to cams. 
Stiffness of the springs and the shape of the cams can be altered to prescribe a wide 
variety of linear and nonlinear structural responses. 

Procedures 

The pitch acceleration was measured by a linear accelerometer, which measured 
accelerations along one axis. The accelerometer was mounted 0.157 m from the rotational 
axis and orthogonal to the y-direction (forward-aft) when the airfoil was at a zero angle of 
attack, giving no acceleration in the y-direction. However, a small portion of the plunge 
acceleration was detected when the airfoil was deflected. The elastic axis was three- 
tenths of the chord length forward of the midchord. 

The free-stream velocity was increased in increments of 2 m/s from 4 to 22 m/s. 
Aeroelastic responses were recorded for 45 s while the free-stream velocity was held 
constant. Flutter was observed to be induced at about 13.5 m/s. Pitch and plunge 
displacements and accelerations of this aeroelastic system were sampled at 525 Hz. 

After recording, the experimental data was decimated by a factor of 2, resulting in a 
final sampling rate of 262.5 Hz. The system was identified using the NARMAX approach, 
as outlined in the section titled “Analysis of NARMAX Model Parameters,” except that 
A/ = 15,600 points was used for estimation and A/^7800 points was used for 
validation. 
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Results 


The results of identifying 10 trials of wing section experiments are presented. 
Figures 5(a) and 5(b) show a typical pitch-plunge displacement and velocity trial used for 
this analysis. The data represents pitch and plunge displacement and velocity sequences 
while the free-stream velocity in the wind tunnel was held constant at 16 m/s. The 
characteristics of this trial are consistent with those reported in previous work (ref. 33). The 
lower panels of figures 5(c) and 5(d) display a 5-s slice of the cross-validation (predicted) 
outputs superimposed on top of the measured outputs for this trial. The predicted outputs 
matched the measured outputs with over 98 %QF. 

Figure 6 shows the cross-validation %QF for each trial. The results show that the 
predicted outputs, for these parameter estimates, account for a large portion of the 
variance. For pitch velocity, the range of %QF is from a minimum of 99.73 percent to a 
maximum of 99.98 percent. For plunge velocity, the range of %QF is from a minimum of 
98.09 percent to a maximum of 99.83 percent. From the 10 trials examined for this study, 
80 percent of predicted outputs accounted for more than 99 %QF of the measured output 
for both pitch and plunge velocity. This result indicates that the NARMAX parameters 
explain the measured data well. Moreover, for every data set the standard deviation 
(STD) of each model parameter was computed at the 95 percent confidence level. These 
results showed that the STDs did not contain zero and suggest that the estimated models 
are accurate. A model parameter whose STD encompasses zero may indicate a spurious 
model term and, hence, should be reformulated (ref. 34). 
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(a). Upper panel: Pitch displacement, (b). Upper panel: Plunge displacement. 
Lower panel: Measured pitch velocity. Lower panel: Measured plunge velocity. 
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(c). Upper panel: Pitch displacement. Lower (d). Upper panel: Plunge displacement, 
panel: Measured and cross-validated pitch Lower panel: Measured and cross-validated 
velocity, 99.98 %QF. plunge velocity, 98.95 %QF. 

Figure 5. (a) and (b). Upper and lower panels: Typical recorded pitch and plunge 
displacement and velocity; (c) and (d). Cross-validation: upper panels: Five-second slice 
of pitch and plunge displacement; lower panels: Five-second slice of predicted pitch and 
plunge velocity outputs superimposed on top of measured velocity output for identified 
NARMAX pitch-plunge velocity models for experimental data set with N =7800. 
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Figure 6. Cross-validation. %QF versus tree-stream velocity (experimental trial). 
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DISCUSSION 

NARMAX Representation of Pitch-Plunge Velocity Dynamics 

The theoretical results demonstrate that the nonlinear difference equation description 
for the pitch-plunge models are NARMAX models. Simulation results show that the 
NARMAX models match the continuous-time response well. These results suggest that 
parametric nonlinear model forms such as the NARMAX class can be used for modeling 
aerospace systems. 

Nonlinear models have the advantage of covering a wider range of system dynamics 
than linear models, which could allow for slower envelope expansion. Using nonlinear 
models to characterize aeroelastic phenomena can provide significant time and cost 
savings for test and development of aerospace vehicles. Moreover, the discrete nonlinear 
models of pitch-plunge provide excellent predictions that could be used for control 
synthesis and statistical studies of NARMAX coefficients and may be of direct relevance 
for health monitoring of aerostructures. 

Discrete-Time Parameter Estimation of Simulated Aeroelastic System 

Simulation studies previously described in the section titled “Analysis of NARMAX 
Model Parameters” showed that, fora NARMAX model representation, the mean of Monte 
Carlo estimates for NARMAX parameters matched the theoretical values well for all SNR 


levels. However, estimates of some parameters (e.g., y 10 ) did not correspond well to 
theoretically computed values. As previously stated, lagged-error terms are difficult to 
identify accurately. Error terms represent a stochastic process that cannot be measured. 
This stochastic process is approximated by a deterministic signal of prediction errors that 
is only an (poor) estimate of the noise (refs. 22 and 28). 

Identification of Experimental Aeroelastic Data 

High %QF cross-validation fits obtained for parameter estimates using NARMAX 
identification methods (fig. 6) show that the identified parameters explain the experimental 
data well. Using %QF alone as an indicator of model goodness may lead to incorrect 
interpretations of model validity. However, in many cases for nonlinear models, %QF may 
be the only indicator that is readily available. 

A model validation technique for nonlinear systems, using higher order correlations, 
was developed by Billings and Voon (refs. 35 and 36). Korenberg and Hunter (ref. 37) 
showed that this model validation technique fails for simple cases. Therefore, this 
approach was not implemented in favor of using the %QF alone as an indicator of model 
goodness. 

When studying aeroelastic systems, assuming that the exact model order and 
structure are well-known a priori may not be practical. In aerospace systems analysis, 
one of the main objectives is not only to estimate system parameters but also to gain 
insight into the structure of the underlying system. Investigation in a future study would 
be worthwhile to determine the result of allowing NARMAX structure detection methods 
(refs. 34 and 38-40) to analyze the data to find the best structure from the data set. 
This may then indicate deficiencies in the analytical model and could lead to improved 
modeling strategies. 


CONCLUSIONS 

Theoretical results demonstrate that the nonlinear difference equation description for 
the pitch-plunge model is a NARMAX model. Simulation results show that the NARMAX 
model matches the continuous-time response well. Moreover, this report contributes to the 
understanding of the use of parametric identification techniques for modeling of nonlinear 
aerospace systems. The main point here is that the NARMAX form is clearly amenable to 
the study of a wide range of aerospace systems and could be computationally efficient. 
NARMAX modeling and identification techniques should be examined further especially 
in the case of severe nonlinear behavior. 
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